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ABSTRACT 

The variability of young stellar objects (YSO) changes their brightness and color preventing a proper classification in traditional 
color-color and color magnitude diagrams. We have explored the feasibility of the fiux variation gradient (FVG) method for YSOs, 
using H and K band monitoring data of the star forming region RCW 38 obtained at the University Observatory Bochum in Chile. 
Simultaneous multi-epoch fiux measurements follow a linear relation = a + j3 • Fx for almost all YSOs with large variability 
amplitude. The slope gives the mean HK color temperature T^ar of the varying component. Because T^ar is hotter than the dust 
sublimation temperature, we have tentatively assigned it to stellar variations. If the gradient does not meet the origin of the flux- 
flux diagram, an additional non- or less-varying component may be required. If the variability amplitude is larger at the shorter 
wavelength, e.g. cr < 0, this component is cooler than the star (e.g. a circumstellar disk); vice versa, if cr > 0, the component is hotter 
like a scattering halo or even a companion star. We here present examples of two YSOs, where the HK FVG implies the presence 
of a circumstellar disk; this finding is consistent with additional data at J and L. One YSO shows a clear ^-band excess in the JHK 
color-color diagram, while the significance of a ^^T-excess in the other YSO depends on the measurement epoch. Disentangling the 
contributions of star and disk it turns out that the two YSOs have huge variability amplitudes (~ 3 - 5 mag). The HK FVG analysis is 
a powerful complementary tool to analyze the varying components of YSOs and worth further exploration of monitoring data at other 
wavelengths. 
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1. Introduction 

The paradigm of a low-mass young stellar object (YSO) in¬ 
volves a pre-main-sequence star surrounded by a circumstel¬ 
lar disk (CSD), a (bi)-polar reflection nebula and an envelope. 
Establishing the presence of a CSD is an observational chal¬ 
lenge. 

The classical technique is to reveal the presence of (T ~ 
1500 K) dust of the CSD by means of a ^-excess in JHK color- 
color diagrams; the YSO is then located at large H - K val¬ 
ues compar ed to small J - H, i.e. right-hand of the reddening 
vector (e.g. |Meyer et~aL]|1997| ). One problem of this technique 
is to distinguish between the various emission components viz. 
star, disk, and a possible reflection nebula. When the CSD is 
seen nearly edge-on, the star is dimmed and heavily reddened 
so that the apparent stellar temperature approaches that of the 
CSD dust. Likewise only a minor contribution from the less red¬ 
dened reflection nebula may be recognisable. In the JHK dia¬ 
gram such a type-2 YSO will be located at both large J - H 
and large H - K values, hence only barely distinguishable from 
highly reddened sources without CSDs (e.g. [Kenyon et al.|1993 
[Stark et al.|2006| ). 

Instead of analysing multi-component YSOs in magnitude 
units, it appears more promising to consider flux (density) or en¬ 
ergy units. Meanwhile numerous multi-wavelength spectral en¬ 
ergy distributions (SEDs) are available, but mostly observed at 
different epochs; this limits potentially the results of SED anal¬ 
yses in case of variable objects. However, particularly YSOs un¬ 
dergo strong brightness variations at X-ray, optical and near in¬ 
frared (NIR) wavelengths, for instance due to stochastic fluctua¬ 


tions of the accretion or to rotating hot/cool spots on the stellar 
surface. 

Therefore a valuable complement to multi-wavelength SED 
studies is the variability monitoring of YSOs, preferably in the 
NIR to be less sensitive to extinction. Simultaneous data in two 
or more Alters provide valuable color information. Assuming 
that the variations originate in the star and not in the circum¬ 
stellar matter, then - in a given aperture - one measures the su¬ 
perposition of a varying hot stellar component and non-varying 
cool disk component. 


This is exactly the analog of what the two-fllter flux variation 
gradient (EVG) method reveals in Seyfert galaxies: a hot vary¬ 
ing nucleus and a cool non-varying contamination by the host 


galaxy ([Choloniewski 198 1[ 

Winkler et al.[[1992[ [Glass 

2004 

[Sakata et al.[[2010[ [Haas et al.[[2011[ [Pozo Nunez et al. 

2012 


2013, 2014, 2015). If the central source is sufficiently hot to al¬ 
low for a Rayleigh-Jeans approximation, the flux is proportional 
to the temperature, and therefore a linear relation between fluxes 
from the two Alters is expected. 

To separate the contributions from variable hot stars and a 
more constant circumstellar environment, we have explored the 
capability of the EVG method for YSOs, using H and K band 
monitoring data of the star forming region RCW 38 obtained at 
the University Observatory Bochum in Chile. Here, we report on 
a first feasibility study, with the focus on the linear relationship 
between the H and K fluxes, and subsequent potential applica¬ 
tions. 
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Fig. 1. Color-color diagram of variable sources in RCW 38. 
Black dots with error bars denote averaged photometric sin¬ 
gle observations from Dorr et al (2013). The main sequence is 
marked in blue and the locus of unreddened classical TTauri 
stars (CTTS) by the red long-dashed line. Two sources with 
ID 557 and 559 are highlighted in red; the small circles mark 
the individual H and K measurements placed at mean J -H due 
to the lack of J variability data; the large circles give the mean 
H - K. The Ay vectors are based on Rieke & Lebofsky| ( |1985| ). 


2. Near-infrared data 


The observations, data red uction and PSF pho tometry have been 
described and tabulated by [Dorr et al. ( |2Q13| ). In bri ef, RCW 38 
ha s been monitored with the 80 cm IRIS telescope ( |Hodapp et| 
al. [2010 ) for three months in spring 2011 with a median daily 


cadence in the /, H, Kg bands (henceforth for short Kg is denoted 
as K). The H and K band light curves are of good photometric 
quality; however, the bright nebula causes a too large uncertainty 
of the J band light curves preventing a robust J band variability 
analysis. 

To analyse the mean SED of each YSO, we used the photom¬ 
etry extracted from the coadded JHK images, reaching about 
1 - 2 mag deeper than 2MASS, and complemented the SEDs 
with Spitzer/IRAC photometry at 3.6 fim from [Winston et al. 

Eig.[^ shows the JHK color-color diagram of 122 variable 
sources with K < 15 mag and having H and K errors smaller 
than 0.1 mag (Table 1 of Dorr et al. 2013). Eor ah but two stars 
the averaged photometric values are plotted. Two stars (ID 557 = 
2MASS J08591359-4733087 and ID 559 = 2MASS J08590708- 


4731499) are selected for a dedicated EVG analysis here. The 
small red circles mark their individual H and K measurements. 
They demonstrate how far a variable source may spread along 
the H - K axis. If J variability data were available, one ex¬ 
pects them to increase the spread of the individual measure¬ 
ments across the J - H ys H - K plane. This demands caution 
when classifying YSOs with single epoch data. While ID 559 
shows a strong K band excess at ah epochs, the evidence for a 
circumstellar disk around ID 557 depends on the measurement 
epoch. If observed only at an epoch of small H - K, ID 557 
could be misidentihed as a main-sequence MO star reddened by 
Ay = 12 mag, instead of a CTTS with mean Ay = 9 mag. 

Our sources ID 557 and ID 559 correspond to Winston et al.’s 
sources number 363 and 325, respectively. Eor both sources, the 


ID=557, J/H/K/L=16.8/15.0/13.8/11.9 mag 



ID=559, J/H/K/L=18.7/16.4/14.4/12.3 mag 



Kflux [mJy] 


Fig. 2. Elux - flux diagram of quasi-simultaneous H ys. K mea¬ 
surements taken with a time separation of less than two hours 
(black dots with error bars). The data are linearly fit by the blue 
solid line, with an uncertainty range indicated by the dotted blue 
lines. The fit parameters are labeled in the upper left corner. The 
red long-dashed line and the red circle mark the constraints for 
an additional non-varying circumstellar disk, derived from the 
color temperature Tyar of the varying star and the mean SEDs 
shown in Eig. The blue dash-dotted lines mark the average H 
and K flux of the star, given by the difference between the black 
and red circle. 


inclusion of the Spitzer photometry argues in favour of an in¬ 
frared excess due to circumstellar dust; Winston et al. assigned 
’’Class 11” to ID 557 and ”Elat-Spectrum” to ID 559. In addition, 
the HK light curves are chaotic with jumps between days sug¬ 
gesting that irregular accretion causes the brightness variations 
via hot spots, and this in turn argues for the presence of a cir- 
cumstehar disk supplying the accretion material. 


3. HIK Flux Variation Gradients 

Eig. shows the flux - flux diagram for the two sources ID 557 
and ID 559. In this diagram, for each (nearly) simultaneous 
HK observation Fh is plotted versus Fk (black dots with er¬ 
ror bars). The variability data clearly follow a linear relation. 
The blue solid line marks the slope p of the fitted linear rela¬ 
tion Fh = a + p ' Fk\ this is the flux variation gradient (EVG). 
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The blue dotted lines delineate the fit uncertainty of the gradi¬ 
ent. The scatter of the data around the gradient is larger than 
expected from the error bars. An explanation could be that the H 
and K observations were obtained with a lag of about one hour 
- hence not sufficiently simultaneous given the strong variations 
from one day to the next. 

The existence of a linear HK flux relationship shows that the 
variable component is intrinsically (i.e. before reddening) hot 
enough to be approximated by the Rayleigh-Jeans tail of a black 
body. The FVG slope gives to first order the HK color temper¬ 
ature, Tyar, of the variable component. The apparent T^ar values 
of 1352 K and 2089 K (Fig. are obviously too low for a stel¬ 
lar source. However, they turn into physically meaningful stellar 
temperature ranges when applying proper dereddening (see Sect. 
1^. Then T^ar hes well above the expected sublimation temper¬ 
ature of dust < 2000K), rejecting the possibility that the 
variations are dominated by temperature changes of the dust. 

To facilitate the illustration of the FVG method, we here ten¬ 
tatively assume that the varying component is the star itself pep¬ 
pered with hot spots (r ~8000K); variations of the CSD, the 
reflection nebula or the envelope are neglected. In Sect. we 
discuss alternatives to these assumptions, however, they do not 
affect the basic conclusions from the FVG technique. 

A striking feature of both YSOs in Fig. is that the FVG 
hits the x-axis at ^ 0.5 mJy and does not pass through the ori¬ 
gin of the diagrams. Note that this is equivalent to the fact that 
the shorter wavelength (F/^) varies with larger amplitude than the 
longer wavelength (^) [^Because the flux variations follow a lin¬ 
ear gradient over a large range, it is unlikely that the gradient 
will change to hit the origin, if the star fades below its currently 
measured H band minimum. Instead, the flux zero point of the 
star might lie somewhere on the gradient, between the lowest 
measured H flux and H = This suggests the presence of an 
additional component which is cooler than the varying star. 

We here identify this additional component with the CSD; 
the envelope is expected to be too cold to contribute at JHK. 
Note that for ID 559 the K band excess already implies the pres¬ 
ence of such a CSD component; for ID 557 the HKL data sug¬ 
gest the presence of warm dust. The location of the CSD in Fig.|^ 
is marked with a red circle, and we justify the CSD color esti¬ 
mates in Section]^ Any CSD variations are either negligible or 
on an uncorrelated time scale different from that of the star, as 
indicated by the small scatter of the data points around the gra¬ 
dient. 

Not only the two selected sources but almost all sources in 
RCW 38 with large variability amplitudes A show a linear re¬ 
lationship between Fh and Fk- About 50% of the 122 variable 
sources have A > 0.4. About 90% of them clearly show a linear 
gradient, which for ~80% of the sources does not go through the 
origin of the flux-flux diagram. 

4, Decomposition of the mean JHKL SEDs 

To constrain in Fig. the zero point of the star on the FVG and 
the contribution of the inferred CSD, we make use of the mean 
SED in JHK and L (i.e. Spitzer/IRAC 3.6 //m). While the L 
band data were measured at a different epoch, one may expect 
that the CSD outshines any stellar contribution at L, hence that 
any variability at L is small compared to JHK. 

Fig. shows the mean SEDs at JHKL. The HK color tem¬ 
perature T^ar of the star as determined from the FVG analysis is 

^ The amplitudes A = {max - min)/avg are A(H) = 0.49, A(^) = 0.32 
for ID 557 and A(H) = 0.84, AiK) = 0.44 for ID 559. 


1D= 557, mean SED 



wavelength [ pm ] 

1D= 559, mean SED 



wavelength [ pm ] 

Fig. 3. Mean Spectral Energy Distributions of ID 557 and 559. 
The total SED (black) is decomposed into a star (blue) and a cir¬ 
cumstellar disk (red). The color temperature of the star is identi¬ 
fied with Tyar as determined from the FVG analysis (Fig.|^. The 
star’s maximum possible mean flux is constrained by the lowest 
total SED, here the J band flux. The remaining SED (total minus 
star, depicted as red diamonds and connected by red solid lines) 
is fitted by a black body (red dotted line). The sum of the star 
and the CSD black body yields the ’’modeled” total SED (black 
dotted lines). 


bluer than that of the total source. For simplicity we model the 
star’s JHKL SED as a black body with temperature T^ar- This al¬ 
lows us to constrain the maximum possible stellar flux, such that 
it matches the total SED at the J band; in this case we assume 
that the J band emission of the CSD is negligible. 

To check whether the remaining SED (total SED minus 
’’maximum possible” star) has a reasonable shape consistent 
with expectations for a CSD, we fit it by a black body leaving 
temperature and intensity as free parameters. As shown in Fig[^ 
the fit (dotted red line) deviates from the data (red diamonds and 
solid lines) by less than 20%. Furthermore, the sum of the star 
and the dust black bodies yields the ’’modeled” total SED (black 
dotted lines). Again, the fit at HKL is good; at J the modeled 
total SED lies slightly above the data, but this is expected from 
setting the star’s flux equal to the total J flux. 

The resulting mean HK fluxes for the CSDs are plotted as 
red circles in Fig. and the maximum possible average star 
fluxes are marked with the dashed-dotted blue lines (difference 
between mean star flux and mean CSD). The true contribution 
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of the CSD and the star may be a bit larger/smaller, respectively, 
because in the above calculation we have assumed that the CSD 
does not emit at J and that the stellar SED follows a black body. 

Of course, there may be alternatives to provide constraints 
on decomposing the SEDs and to fit a star and a CSD. The in¬ 
vestigation of alternatives is postponed to the future. 


5. Dereddening of star and disk 

In Eig.j^the derived HK color temperature of star and CSD are 
about Tstar ~ 1300,2100K and Tdust ~ 800, HOOK. This is 
much lower than the range expected for an unreddened PMS star, 
with Tstar ~ 3000 - 8000 K and hot dust {T^ust ~ 1500 K) in the 
CSD. Thus, we check if reddening has ’’shifted” Tstar and Tdust 
from reasonable temperatures to the observed low values. 

We assume that both star and CSD are reddened by the same 
amount of extinction. Dereddening the JHKL photometry by the 
Ay values from Pig. and the extinction curve from |Rieke &| 
[Lebofsky ( |1985| ) we obtain dereddened fiux-fiux diagrams that 
look qualitatively similar to those in Pig.|^but where the scaling 
differs in the sense that Fh is more stretched than Fk- Por Ay = 9 
the dereddened values for ID 557 turn into Tstar ~ 5800 K and 
Tdust ~ 1350 K; an Ay value of 10 mag yields Tstar ~ 7700 K 
andr^,,, ^ 1400K. 

Notably, dereddening of ID 559 with Ay = 10 (Pig.[2) yields 
Tstar ~ 2300 and Tdust ~ HOOK, hence does not sufficiently 
raise the star and dust temperatures. However, increasing the 
visual extinction to Ay = 17 mag yields Tstar ~ 4900 K and 
Tdust ~ 1350 K. In the JHK color-color diagram dereddening 
with Ay = 17 shifts ID 559 below the CTTS line (Pig.[^. This 
’’apparent” inconsistency could be explained by an unresolved 
reflection nebula diluting the actual reddening (e.g. [Stark et al. 
20061 |Krugel]200^ . 


6. Discussion 

Given the complexity of YSOs, the different components may 
contribute in a different manner to the variability observed across 
the optical, NIR and MIR wavelength range. Based on a simul¬ 
taneous optical {R ba nd) and MIR (IRAC 3.6-4.5 jim) variability 
study of NGC 2264, [Cody et ah] ( |2014| ) present an interesting 
discussion on different possible causes of the variability. Por es¬ 
sentially all disk-bearing sources the variations at both IRAC fil¬ 
ters are correlated. These variations are also correlated with the 
R band variations in 40% of the sources, indicating increased 
heating of the disk in response to variable accretion or hot spots 
on the surface of the star. Por the remaining 60% of Cody et 
al.’s sample with uncorrelated optical and MIR variability, we 
propose an alternative explanation. Young stars are believed to 
be bom as binaries or in multiple systems. If of different mass, 
the two binary components evolve differently, resulting in a blue 
star and an (infra)-red star. They may show independent varia¬ 
tions, the blue star dominating at optical and and the red star at 
MIR wavelengths. If such a binary is unresolved, an object with 
uncorrelated optical and MIR variability will be observed. 

Swirling dust clouds (moving across the line of sight towards 
the star) produce varying extinction, resulting i n a larger vari- 
ability amplitude at the shorter wavelengths (e.g .|Skmtskie et al!] 
|1996[ [Carpenter et al.|[200T| [Indebetouw et al.[[2006[ ). Varying 

extinction may explain the variability behaviour of UX Orionis 
stars believed to be seen at grazing angles of the disk (e.g. [Herb^ 
et al.|1994^[Bertout|2000[[Dullemond et al.|2003[ ). If varying ex¬ 


tinction is the dominant mechanism, then in a JHK color-color 


diagram or a color-magnitude diagram (CMD) the varying data 
points of each object should be aligned with the extinction vec¬ 
tor. Such an alignment is only rarely found, e.g. in less than 10% 
of both the Ori on data of Carpenter et al. and our M 17 data 
( [Scheyda |2010[ ). Por RCW 38, light curves are available only at 
H and K, but not at J. Our two YSOs (ID 557 and ID 559) show a 
clear elongation in the two CMDs HI{H-K) and K/(i7-K), which 
however is not well aligned with the direction of the extinction 
vector|^We suggest that for both sources variable extinction may 
be present but plays a minor role. 

As shown in Pig.[^ the CSD contributes by about 50% to the 
total H and K flux; this holds also for the dereddened data. Thus, 
after subtraction of the CSD, the two young stars have a huge 
variability amplitude (A ~ 20 for ID 557 and ID 559 reaches 
A ~ 100, equivalent to ~5 mag). This indicates that their lumi¬ 
nosity is essentially powered by accretion bursts. Puture studies 
may provide clues to whether such strong flux variations can be 
explained only by temperature and/or area changes of hot spots, 
or whether there are contributions by a variable hot gaseous disk, 
in addition to the cooler circumstellar disk or whether the huge 
amplitudes are caused by other effects. 

So far, with the aim to facilitate the illustration of the PVG 
method, we have neglected variations of the CSD, the reflection 
nebula or the envelope ^Now we address such variations. If there 
is variable heating of the dust disk in response to variable accre¬ 
tion or hot spots on the surface of the star, and if the dust reacts 
fast enough (e.g. simultaneous within the duration of the H and 
K observations), then the PVG measures a combination of the 
hot stellar and the cool dust flux variations. Then (after dered¬ 
dening) T^ar is not the temperature of the varying stellar compo¬ 
nent ~ 8000 K), rather it is lowered by the contribution 

of the varying dust component (T^^f ~ 1500K); hence Tyar un¬ 
derestimates Accounting for the constant star and dust 

contributions and using T^ar instead of in the decomposi¬ 

tion of the mean SED (Pig.[^ may result in an underestimation 
of the stellar flux at H and K, and hence an overestimation of the 
variability amplitude. A comprehensive treatment of the many 
combinations and aspects of NIR PVGs for YSOs, including the 
effects of reflection nebula and envelope, is planed for the future. 


7. Summary and outlook 


We have explored to which extent the flux variation gradient 
(PVG) technique can be transferred from active galactic nuclei 
to varying young stellar objects (YSOs). The results from HK 
monitoring data of the star forming region RCW 38 are: 

• Por almost all YSOs with large variability amplitudes the 
flux variations follow a linear relation Fh = cc+P'Fk- Such 
a linear relation between fluxes from the two filters is con¬ 
sistent with the Rayleigh-Jeans approximation. In this case 
the PVG slope p gives the color temperature of the varying 
component. 


^ The resulting ratios AhIAk are 1.84 + 0.09 and 2.19 + 0.10 for 
ID 559 and ID 557, respectively, which is 3.1 cr and 6.2 cr apart from the 
standard extinctio n law {AhIAk = 1.56, with a ratio of total-to-selective 
exti nction R - 3.1,[Rieke & Lebofsky[1985| ). A much larger value ~ 4 
(e.g. [Hoffmeister et al.[2008| ) would be reqired (locally at the sources) to 
align the H and K variability of ID 559 and ID 557 with the direction of 
the extinction vector. Such a large R value does not fit the overall con¬ 
sistency of the RCW 38 sources in the JHK color-color diagram with 
the standard extinction law (e.g. Fig. 3 in Dorr et al. 2013|). 

^ The assumption that the variations occur also in the CSD, implies 
the existence of a CSD. Hence, in this case the FVG technique reveals 
the presence of a CSD. 


4 









































R Pozo Nunez et al.: Circumstellar Disks and the FVG method 


• Temperature considerations indicate that the star dominates 
the flux variations and that variations of the circumstellar 
disk/envelope or the reflection nebula play a minor role. 

• Two examples where the variability amplitudes are larger 
at shorter wavelenghts and where variable extinction might 
play a minor role, are investigated in more detail. The neg¬ 
ative FVG offsets (a < 0) imply the presence of an addi¬ 
tional cool component such as a circumstellar disk. At least 
for some cases the FVG technique appears to be more sen¬ 
sitive for detecting circumstellar disks than the ^-excess in 
the JHK color color diagram. It appears worth using FVGs 
to reexamine the frequency of disks and the evolution of the 
disk dispersal in star forming clusters. 

• The low apparent HK color temperatures for the star and the 
disk, as inferred from FVGs, can be explained by substan¬ 
tial reddening. The extinction, required to deredden star and 
disk to realistic intrinsic temperatures, is at least as high as 
inferred from JHK color color diagrams, consistent with ex¬ 
tinction and scattering models. 

• FVGs provide basic input for the decomposition of SEDs. 
A more stringent SED decomposition requires further con¬ 
straints which may be obtained from additional monitoring 
in a third Alter and subsequent EVG analysis across several 
Alter pairs. 

• After subtraction of the circumstellar disk, the two young 
stars investigated here display a huge variability amplitude 
(~ 3 - 5 mag). If the high amplitudes are not pretended by 
other effects to be explored, the luminosity of these YSOs is 
essentially powered by accretion bursts. 

To conclude, EVGs are a powerful tool not only to infer the pres¬ 
ence of a circumstellar disk but also to disentangle stellar and 
disk parameters with the help of reasonable assumptions. 

To give an outlook on further applications, we have tenta¬ 
tively applied the EVG technique to two data sets: 






Simultaneous JHK monitoring of the star forming region 
M 17 with the Infrared Service Eacility (IRSE) at the 1.4 m 
telescope in Sutherland, South Africa ( Scheyda |2010| ) yield 
for both Alter combinations, JH and HK, that the flux vari¬ 
ations follow extremely well a linear relation; Tyar agrees 
within the errors. Note that JK is redundant, being the prod¬ 
uct of JH and HK. The low scatter around the EVG is prob¬ 
ably also a consequence of the perfect simultaneity of the 
observations. There are YSOs where the HK EVG with neg¬ 
ative a implies a circumstellar disk; simultaneously the JH 
EVG has a positive a > 0 implying a reflection nebula or - 
more consistent with the unresolved appearance on optical 
BVRl images - an unresolved blue companion star. Since 
YSOs are suspected to be born in multiple systems, such a 
companion is not unexpected, but needs to be verified. Eor 
instance, the young PMS binary XZ Tau consists of a rather 
evolved blue star and a less evolved infrared star (|Haas et al.| 
p[990l ). 


The YSOVAR project lists Spitzer/IRAC 3.6 a nd 4.5//m 
monitoring photometry of IC 1396 A and Orio n ([Morales- 
[Calderon et al.||200^ [Morales-Calderon et al.||2011| ). 387 
Orion sources have amplitudes A > 0.2, and we find that 
their 3.6 and 4.5 jjm flux variations follow a linear relation. 
Eor 154 of these sources (40%) we find Tyar > 1500 K, and 
for 64 sources (16%) Tyar > 2000 K, reaching 3500 K. 

Given that the sources are reddened, the variable component 
of these sources appears too hot for being dust. On the other 
hand, at 3.6 and 4.5 /im most of the e mission is due to the in¬ 
ner disk and envelope, if present (e.g. IRebull et al.|20T^ .To 


bring both findings into a consistent picture, both stars and 
disk/envelope may vary simultaneously with similar strength 
so that the EVG measures an average Tyar, higher than for 
dust but slightly lower than for the star. 

These two data sets {JHK of M 17, and YSOVAR) strengthen 
the suggestion that the EVG technique opens the door to valuable 
insights into YSOs which are hard to obtain otherwise. 
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